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We present a family of nonautonomous bright and dark soliton solutions of Bose-Einstein conden- 
sates with the time-dependent scattering length in an expulsive parabolic potential. These solutions 
show that the amplitude, width, and velocity of soliton can be manipulated by adjusting the atomic 
scattering length via Feshbach resonance. For the cases of both attractive and repulsive interaction, 
the total particles number is a conservation quantity, but the peak (dip) density can be controlled 
by the Feshbach resonance parameter. Especially, we investigate the modulation instability process 
in uniform Bose-Einstein condensates with attractive interaction and nonvanishing background, and 
clarify that the procedure of pattern formation is in fact the superposition of the perturbed dark 
and bright solitary wave. At last, we give the analytical expressions of nonautonomous dark one- 
and two-soliton solutions for repulsive interaction, and investigate their properties analytically. 

PACS numbers: 03.75.Lm, 05.30.Jp, 67.40.Fd 

Keywords: Nonautonomous soliton solution; interaction; Bose-Einstein condensation 



I. INTRODUCTION 

The classical soliton concept was introduced firstly by 
Zabusky and Kruskal [l| for autonomous nonlinear and 
dispersive dynamic systems where the time variable has 
only played the role of an independent variable and has 
not appeared explicitly in the coefficients of the non- 
linear evolution equation. These autonomous soiitons 
do not disperse and completely preserve their localized 
form and speeds during propagation which has motived 
a great attention in optical fibers and condensate physics. 
As a theoretical model, the nonlinear Schrodinger equa- 
tion has been adopted extensively to govern the dynamics 
of autonomous bright and dark soiitons in optical fibers 
0, [|[ and Bose-Einstein condensates (BECs) [ij . 

In BECs bright and dark soiitons have been paid more 
particular interest experimentally and theoretically. The 
bright soliton [Hj-FTTj is expected by the balance between 
dispersion and attractive mean-field energy. However, 
collapse of bright soliton [l2[ may occur owing the at- 
tractive interaction of bosons. To avoid this collapse, one 
should restrict the BECs dynamics into the quasi-one- 
dimensional regime, i.e., the energy of two body interac- 
tion is much less than the kinetic energy in the transverse 
direction. The dark soliton [l3l - [l7j can be formed in the 
case of repulsive interaction of bosons. It denotes the 
macroscopic excitation characterized by a local density 
minimum, and a phase gradient of the wave function at 
the position of the minimum. Comparing with the at- 
tractive case, large condensates can be realized for dark 
soliton owing the repulsive interaction of bosons. 

When the physical system is subjected to various exter- 
nal time-dependent forces, the nonautonomous nonlinear 
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evolution models typically arise and the term of nonau- 
tonomous soiitons was introduced firstly. In fact, 
different aspects of dynamics in nonautonomous mod- 
els [l9l - [2l| and the controllable soliton solutions [HI, [22| — 
HU in optical fibers have been investigated theoretically. 
For BECs the nonlinearity resulting from the interatomic 
interaction is denoted by the effective scattering length 
which can be tuned experimentally by utilizing the Fes- 
hbach resonance [25| . even including its sign. A sinu- 
soidal variation of the scattering length has also been 
used to form patterns such as Faraday waves [26|, [27j . 
The controlling soliton trains [28| was also created start- 
ing from periodic waves. Therefore, the BECs soiitons 
forming by magnetically tuning the interatomic interac- 
tion near a Feshbach resonance are a typical example of 
such nonautonomous soiitons in external potentials and 
offer a good opportunity for the nonlinear excitations ex- 
ploration [29j. The properties of such nonautonomous 
soiitons in BECs are not well explored and it is our pur- 
pose in the present paper. 

At the mean-field level, the evolution of the macro- 
scopic wave function of BECs can be described by the 
Gross-Pitaevskii equation [3(J HH , 

F)^i) (r t\ h 2 v 2 

where $ (r,£) is normalized to the number of condensed 
atoms, i.e., N = J |\?| d 3 r, m is the atomic mass, V cx t 
is a harmonic trap given by V ext (r) = m/2[uj 2 l x + 
w ± (y 2 + z 2 j] with wo and ui± being the axial and trans- 
verse harmonic oscillator frequencies, and the effective 
interatomic interaction reads g = ^■nh 2 a s /m with a s 
being the s-wave scattering length (a s < for attrac- 
tive interaction; while a s > for repulsive interac- 
tion). The linear oscillator lengths in the transverse and 
cigar-axis directions is defined by l± = (h/muj±) 1 ^ 2 and 
lo = (h/m l^ol) 1 ^ 2 ! respectively. For a cigar-shaped con- 
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densate at a relatively low density, when the energy of 
two body interactions is much less than the kinetic en- 
ergy in the transverse direction, i.e. when e 2 = ~ 
N\a s \/l < 1 [IH, where £ = (87m |a s |)~ 1/2 is the 
healing length and n oc N/ (l±lo) is a mean particle 
density, then the system becomes effectively quasi-one- 
dimensional regime [!, HH with the explusive potential 



d 



d 2 



ij^ijj (x, t) = -jj-£il)(x, t) - -X 2 x 2 ip(x, t) 
+ 2N^\iP(x 7 t)\ 2 iP(x,t), 



(2) 



where A = 2 |wo| /wj_ <C f , and the time t and coordinate 
x has been measured in units 2/uj± and l±, respectively. 
In real experiment for BECs Q the scattering length a s 
is tuned with a Feshbach resonance from repulsive to at- 
tractive. For soliton created with the particle number 
N ~ 10 3 , wx = 2?r x 710Hz, u; = 2ni x 20Hz, the pa- 
rameter A = 0.056 is small. With the above conditions, 
the units l± is about lAAfim and the units time for a 
BECs trapped with the transversal size order l± corre- 
sponds to 4.5 x 10 _4 s. The lifetime of a BECs is of the 
order of Is, which is about 220 in our dimensionless units. 
If one chose the scattering length is increased in the form 
of a s (t) = aoexp(Xt) with the initial scattering length 
ao = —0.02 nm. After 40 dimensionless units of time, the 
value of the atomic scattering length turns to a s = —0.19 
nm corresponding to e 2 ~ 0.02 which provides the safe 
range parameters. 

In terms of the transformation ip (x, t) = 
q(X 7 T)exp(\t/2-i\x 2 /4) with X = a; exp (At) 

and T = 2 f* exp (2Ar) dr, Eq. © reduces to the 
standard form 

^ + ^' (Arao/Ml9|2<Z = - (3) 

This result show that in BECs the nonautonomous soli- 
tons formed by magnetically tuning the interatomic in- 
teraction via Feshbach resonance can be obtained from 
the autonomous solitons. In this paper, we will ex- 
plore the generalized nonautonomous bright and dark 
solitons of Eq. As an example, the nonautonomous 
bright soliton solutions are obtained on the vacuum state 
background and nonzero background, respectively. With 
these exact solutions the corresponding dynamic proper- 
ties are discussed in detail. At last, we investigate the 
dynamic behavior of the nonautonomous dark soliton so- 
lutions. 



II. NONAUTONOMOUS BRIGHT SOLITON 
SOLUTIONS 

In this section, we consider the case of attractive inter- 
action between atoms, i.e., ao < 0. In this case, Eq. ([3]) is 
a integrable model, and its soliton solutions can be con- 
structed by several technique, such as inverse scattering 
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3], Darboux transformation [34-36] and Hirota 
37| . With the expulsive parabolic potential and 
the Feshbach-managed time-dependent scattering length, 
we will present the exact nonautonomous bright soliton 
solutions analytically on the vacuum state background 
and nonzero background for Eq. ©. 



A. Nonautonomous bright soliton solution and 
soliton interaction on the vacuum state background 

Firstly, we can present the nonautonomous one-soliton 
solution for Eq. ([2]) as follows 



if) = A s sech# s exp (iip s + At/2) 



(4) 



where 9 S and ip s are given by 



% = fiA s ( xe At - 2k s j e 2AT dr ) - 6 , 



r - = l ,xe xt - ^- + (fi 2 A 2 s - k 2 ) e 2XT dr - (5) 



here /x = yjN \a \ Jl±. The solution tp in Eq. ((H) de- 
scribes a bright soliton of BECs with time-dependent 
atomic scattering length in an expulsive parabolic poten- 
tial, with the initial maximum amplitude A s , the initial 
wave number k s , the initial location 9q/(jj,A s ), and the 
initial phase fo/k s . When A = 0, Eq. ((U) can reduce 
to the solution for the standard nonlinear Schrddinger 
equation. From Eq. ([5]) we see that the soliton in Eq. 
(|4]) can undergo compressing effect with the increasing 
of the scattering length tuned by the Feshbach reso- 
nance. This character does not exist for the situation 
of the uniform nonlinear Schrddinger equation Q . This 
nonautonomous soliton has an increase in the peak value, 
while the particles number is a conservation quantity due 
to /_ °° |^»i_ so i| 2 cfcr = 2A s j[L. The velocity of soliton is 
affected by the Feshbach resonance parameter A, which 
reads V s = -2/c s (e A *-Ae~ At J Q e 2Ar dr). This result show 
that the size of bright soliton in BECs can be tuned by 
adjusting the attractive interactions even in the presence 
of the expulsive potential. 

Another interesting problem is to discuss the interac- 
tion of two nonautonomous bright solitons. To this pur- 
pose we present two-soliton solution of Eq. (|2|) as follows 



exp (At/2) . 



(6) 



where 



F b = fx cosh {9 1 + 6 2 ) + h cosh {0! - 6 2 ) 

+ /3COS (<pi - ip 2 ) , 
G b = gi coshflie^ 2 + g 2 cosh 

+ ig 3 (sinh6»ie^ 2 - sinhfee^ 1 ) , 
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fi = (fc 2 -fcif +M 2 (^2-^if, 
h = (k2-k 1 f+ f i 2 {A 1 +A 2 )\ 
h = -V^M 2 , 

gi = 2A 2 [(k 2 -k l ) 2 +p? (A 2 -A 2 )}, 
,g 2 = 2A 1 [(k 2 - h) 2 + n 2 (Al-At)}, 
g 3 = 4fj,AiA 2 (k 2 - fci) , 



with the parameters 



9j = fj,Aj I xe 



„2At 



dr ) - °0j, 



k 3 xe xt - ^- + (V 2 A 2 - k 2 ) f e 2Xr dr - Voj , (7) 
4 Jo 



here j = 1,2, 0(y and tpoj is an arbitrary real constant, re- 
spectively. The solution in Eq. ([6]) describes a general in- 
teraction between two nonautonomous solitons with the 
different center velocity V s ,i and 14,2, respectively. From 
Eq. (f7]) we get the velocity of each soliton 



V s ,j = -2k sJ (e xt - Xe 



„-A* 



„2Ar 



dr), j = l,2. 



In order to understand the nature of two nonautonomous 
solitons interaction, we analyze the asymptotic behavior 
of two-soliton solution in Eq. ([6]). Asymptotically, the 
solution in Eq. ([6]) can be written as a combination of 
two one solutions in Eq. The asymptotic form of 

two-soliton solution in limits t — s- — oo and t — > oo is 
similar to that in Eq. ((4]). 

(i) Before collision (limit t — > — oo) 

(a) Soliton 1 (0i « 0, 9 2 -oo) 



2-sol 



2^/JT^cosh(0 1 + xo)' 
(b) Soliton 2 ( 2 « 0, 0i -> oo) 



V>2- 



sol 



j 2 e xt / 2 e *(v2-<fe) 



V/1/2 cosh (0 2 - x ) 



(8) 



(9) 



(ii) After collision (limit < — > 00) 
(a) Soliton 1 (0! « 0, 2 -> 00) 



^2- 



sol 



71 e 



At/2 



2V/1/2 cosh (0i - x )' 
(b) Soliton 2 ( 2 « 0, 0i ^ -00) 



2 -sol 



2 V / 7u/2COsh (02+Xo) 



(10) 



(11) 



where 



0i = arctan (33/32) ,02 



arg (g 3 /5i) , 



71 = \/g 2 + sf> 72 = \jg\ +gh x o = ^ In (/2//1) 



From the above asymptotic behavior of two-soliton solu- 
tion, we know that there is no change of the amplitude 
for each soliton during the process of collision. However, 
from Eqs. © to (ITT]) we find there is a phase exchange 
2(f) j,j = 1, 2, and center shift 2xo for soliton 1 and soliton 
2 during collision, respectively. These results show that 
the collision of two nonautonomous solitons is elastic. 



B. Nonautonomous bright soliton solution on 
nonzero background 

It is easy to find two basic solutions of Eq. @. One is 
■0 = 0, which corresponds to the vacuum particle density 
state. The nonautonomous soliton solutions in Eqs. (U]) 
and © are constructed on this zero background to the 
moment. The other interesting solution of Eq. © is a 
plane wave solution 



tp c — ^4exp (iifc + Xt/2) . 



(12) 



with the initial amplitude A, wave number k, and ip c = 
-\x 2 /4: + kxe xt + (2(i 2 A 2 - k 2 ) / Q * e 2XT dr. This solution 
can be seen as the background with the temporal vari- 
ation particle density. It will be very interesting to get 
the exact soliton solution and its dynamic properties on 
such background. Employing Darboux transformation 
[34l - [3"6l | . we obtain the nonautonomous soliton solution 
as follows 



A + A„ 



' cosh + a cos ip 



(13) 



where 



G c = b\ cosh + cos tp + i (b 2 sinh + c sin tp) , 

= M lX e xt - [fxA s M R + (k + k s ) Mj] f e 2Xr dr - O , 

Jo 

Lp = M R xe xt - [{k + k s ) M R - iiA 8 Mi] / e 2Xr dr - tp , 

Jo 



with the parameters b\ 



-2n 2 AA s /D, b 2 



-2fj,AM R /D, c = Mi/(nA„), D = (i 2 A 2 s + M R , and 
Mr + iMj = [(k -k s - ifiA s ) 2 + i^A 2 } 1 / 2 , which imply 
that Mj = as A s = 0. Here 0o, tpo, A s , k s , A, and k is 
an arbitrary real constant, respectively. From Eq. (|13[) 
one can see that, as A vanishes, the solution (fT3")) reduces 
to one soliton solution in Eq. Q. On the other hand, 
when the initial amplitude of the soliton A s vanishes, the 
solution ip in Eq. (|13[) reduces to the solution in Eq. (|12p . 
Therefore, the exact solution (fT3)) describes generally the 
dynamics of the nonautonomous bright soliton embedded 
in the temporal variation particle density background, 
characterized by the envelop propagation velocity 
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Xe~ tx J* e 2Xr dr) 



V sc = {^A.Mr/Mi + k + k : 
Xe~ tx 9o/Mj. 

Based on the above exact solution, we analyze in detail 
the modulation instability process and the formation of 
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spatial pattern in BECs. As discussed in Ref. |34l. l35j. 
we firstly consider a special case, i.e., k s = k. In this 
situation, there are two representative results: 

(i) When 4 A 2 > A 2 , we have Mj = which implies 
that the soliton velocity V sc becomes infinite. It is to 
say that the soliton in Eq. (fT3"|) is completely trapped in 
spatial direction and undergoes the modulation instabil- 
ity process [3(. Indeed, by introducing a small quantity 
e = exp(# ) for 6>o < 0, and then linearizing with respect 
to e we have the approximation for the initial value 

ip(x,Q) « [p + excos^Nhx - <p )]e- tXx2 ' A+lkx , (14) 



where p — (2A 2 
X = A a M\(M\ - 



-A 2 - iA s M 1 )/{2A) with \p\ = \A\, 
iA s )/(2A 2 ^ 



and Mi = y/4A 2 - A 2 . 
The solution of the initial value problem of Eq. @ can 
be well described by the solution (fl3l) under the case of 
Mi = 0. As a result, a small periodic perturbation of the 
plane wave solution may lead to the onset of instability, 
(ii) When AA 2 < A 2 , the soliton (JTSJ) becomes 

*=(-A + M 2 f C0S // Sln ^ e* At+ *% (15) 
\ A s cosh 6 - 2 A cos ip J v ' 

where 9 = pM 2 {xe xt - 2k f e 2Xr dr) - 6 and 

ip = p 2 A s M 2 /o e 2Xr dr - ip with M 2 = y/A% - 4A 2 . 
On the nonzero background the nonautonomous bright 
soliton possesses the same properties as that on zero 
background. The peak value increases and the soliton 
width compresses with the increasing value of the 
scattering length, respectively. The dynamic soliton 
evolution exhibits the periodic oscillation of the am- 
plitude and breather behavior due to the presence 
of nonvanishing background [HI, |36j]. On the other 
hand, under the effect of the expulsive parabolic 
potential, the bright soliton can propagate in the 
longitudinal direction, instead of oscillation in attrac- 
tive parabolic potential. The atoms of the bright 
soliton against the background can be obtained, i.e., 

/-T (lV<| 2 - IV» (±oo, t)?) dx = A s (b 2 + c 2 ) I \Mj\ 

J-°° (A s — 2A cosh 6 cos p) / (cosh 6 + a\ cos p) 2 d6, 
where B = M R /Mj, <p = BO + A, and 
A = fj,A s Mi (B 2 + l) j£ e 2XT dr + B9q — (po- By 
numerically verifying this integration we find it is a 
conservation quantity which indicates that during the 
process of the compression of the bright soliton the 
number of atoms in the bright soliton keeps invari- 
ant. Another interesting problem is how such soliton 
can be created. From the expression (|15[) , we can 
see that the initial wave function can be written as 
the form <f> = -Ae itp " ± iM 2 e itp " sech ([iM 2 x - O ) as 
po = ±7r/2, ±37r/2, • • • . This result shows that the 
solution dTB]) can be generated by coherently adding in 
quadrature a bright soliton to the background. 

To better understand the properties of nonautonomous 
bright soliton on nonvanishing background with the case 
k s k, we can decompose Eq. (|13p into the form 



ip = tp d + ip b , 



(16) 



where 



ip d = [A + A s h + 



ip b = A s e xt ' 2+iv 



iA s b 2 sinh(9 
cosh# + ai cos ip 



\t/2+i Vc 



(l — b 2 ) cos p + ic sin p 
cosh 8 + bi cos ip 



As A = we have ipd = and ipb = A s exp(— iip) sech8 s , 
which shows that in Eq. (1161) the former is zero solution, 
and the latter gives rise to one soliton solution for Eq. 
@. With the increasing of \A\, a dip starts to occur for 
ipd describing a perturbed grey solitary wave under the 
effect of nonzero background. At the same time ipi, ex- 
hibits the periodic oscillation of the amplitude in prop- 
agation which can be considered as a perturbed bright 
solitary wave. Therefore, the solution (IT3l) describes the 
superposition of the perturbed dark and bright solitary 
waves expressing the procedure of the pattern formation. 



III. NONAUTONOMOUS DARK SOLITON 
SOLUTIONS 

In this section, we consider the case of repulsive inter- 
action between atoms, i.e., do > 0. Using the results in 
Ref. [371 . |38| , we obtain the nonautonomous dark soliton 
solution of Eq. ([2]) in the form 



ip = Z-q(x,t)exp(iip d + Xt/2), (17) 
A 4 

where p d = -\\x 2 + £ Q xe xt - (fg + 2n) J* e 2Xr dr + Co 
and q(x, t) is to be determined. For nonautonomous one- 
soliton solution, q(x, t) is given by 

q (x, t) = ~[(1 + Z x ) - (1 - Z x ) tanh |-], (18) 

where 

m = P lX e xt + (^Ak-P 2 - 2£ )Pi j e 2Xr dr + £, 



^Ak - P 2 + iPj. 
- Pi - *Pi ' 



(19) 



From Eq. (|19l) we see that the existence of nonau- 
tonomous dark soliton implies the condition k > P 2 /4. 
From the solution (fl~7| with Eq. (fT8|) we clear two spe- 
cial cases, i.e., Z\ = ±1. When Z\ = 1, the nonau- 
tonomous dark solitonin Eqs. (flTl) reduces to plane 
wave solution ip = ^fR,/ p,e lipd+xt ' 2 , which corresponds 
to the variation distribution density of bosons due to 
the existence of the Feshbach resonance parameter A. 
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When A = 0, this case corresponds to the uniform dis- 
tribution density of bosons. On the other hand, when 
Z x = -1 the solution in Eqs. (17]) and (JTEJ) be- 
comes tpi = — yfR,/ [ie lVd+xt / 2 tanh (771/2), where 771 = 

2y/~Kxe xt — 4i/k£o / * e 2AT <fr + £i. This solution represents 
black soliton solution in BECs with variational amplitude 
and the soliton velocity V d = 2^ (e At -Ae~ At / * e 2(Ar) dr). 

The result in Eqs. (JT7J and ([Tj§ shows that this 
nonautonomous dark soliton undergo compressing ef- 
fect with the increasing of the scattering length, while 
the peak value has an increase. From Eqs. (|17p and 
OH) we have |d min 



= tJk - Pi/4/fJ, exp (Xt/2) which 
is inverse proportion to the initial soliton width 1/P\. 
The dark soliton possesses the accelerated motion with 
the absolute increasing scattering length tuned by the 
Feshbach resonance instead of oscillation in attractive 
parabolic potential and the velocity reads Vd — (2£o — 
>/4k- P 2 )(e At - \e~ xt f* e 2(XT ^ dr) obtained from Eq. 
(UHJ). From Eq. (fTT)) we obtain the particles number in 

the form (\ip\ 2 - \ip (±00, i)| 2 ) dx = -P1//J 2 , which 

is a conservation quantity. The above analysis implies 
that we can control the dip matter density of dark soli- 
ton by adjusting appropriately the Feshbach resonance. 

For nonautonomous dark two-soliton solution, the ex- 
pression of q(x, t) is given by 

. 1 + Z x e* + Z 2 e^ + A 12 Z 1 Z 2 e^+^ 
92 (*»*) = 1 , , ^ , 3 „ n ,4- n , ' ( 2 °) 



where 



jji,- = P 3 a;e At + (yJ^K - P 2 - 2^)P, J e 2Xr dr + Ci,,, 



4k - P 2 + iPj 



4k - P 2 - iPj 



A, 



4k - P X P 2 - y/4n- P 2 ^4k- P 



4k + P 1 P 2 - y/iK- P'I^/Ak- P 



(21) 



solution in Eq. (|T7|) with Eq. (|20l) . Asymptotically, 
the two-soliton waves in Eq. (|2Tj)) can be written as a 
combination of two one-soliton waves in Eq. (ITS)) . The 
asymptotic form of two-soliton solution in limits t — > — 00 
and i — > 00 is similar to that of one-soliton in Eq. (fl7|) 
with Eq. CTj§. 

(i) Before collision (limit t — > —00) 

(a) Soliton 1 (771 w 0, 772 — > —00) 

V-^^y- [1 + £ 1 -(1-Z 1 ) tanh |-], (22) 



(b) Soliton 2 ( ry 2 « 0, 771 — !• 00) 



2// 



[1 + Z 2 - (1-Z 2 ) tanh 



(ii) After collision (limit t —> 00) 
(a) Soliton 1 (771 ~ 0, 772 — >■ 00) 



(23) 



- ^f 72 !! + *-(!- Z0 tanh *±*] 



2/y 

(b) Soliton 2 ( 773 w 0, 771 



(24) 



V^^y- [1 + Z 2 -(1-Z 2 )tanh|-], (25) 

where the center shift of dark soliton is given by So — 
ln^4i2. By analyzing the asymptotic behavior of two- 
soliton solution in detail, we know that there is no change 
of the amplitude for each soliton during collision, while 
one should notice that the factor \Zj\ = 1, j = 1,2, again. 
However, from Eq. (|22|) to Eq. (|25|) we find a phase 
exchange So for soliton 1 and soliton 2 during collision. 
These results show that the collision of two dark solitons 
is elastic. 



IV. CONCLUSIONS 



here j = 1,2. From the solution in Eq. (|20)) , we can see 
its asymptotic behavior 

q 2 (x, t) — > Z X Z 2 , as x — > +00, 
q 2 (x, t) — > 1, as x — > — 00, 

which shows only a phase shift Sd.i + Sd >2 , here Sd.j — 

arctan[2Pj ^4k - P 2 / (4k - 2Pf)],j = 1, 2, as from x -> 

+00 to X -> -00. The solution in Eq. (|T7|| with Eq. (|2TJ|) 
describes a general scattering process of two dark solitary 
waves of BECs on the nonzero background, characterized 
by the different center velocity V\ and V 2 , respectively. 
From Eq. (f2T|) we get the velocity of each soliton as 



Vj = (2£o - jAn-P 2 ){e xt - Xe-^f'e^dr), j = 



JO 

1,2. In order to understand the nature of two solitons 
interaction, we analyze the asymptotic behavior of the 



In summary, we present a family of nonautonomous 
soliton solutions of BECs with the time-dependent in- 
teratomic interaction in an expulsive parabolic potential. 
Our results show that the amplitude, width, and velocity 
of nonautonomous soliton can be affected by the time- 
dependent atomic scattering length, which can be tuned 
by the external filed from the so-called Feshbach reso- 
nance technique. These results also provide an experi- 
mental tool for investigating the range of validity of the 
one-dimensional Gross-Pitaevskii equation with an ab- 
normal expulsive parabolic potential. For the cases of 
both attractive and repulsive interaction, the total parti- 
cles number is a conservation quantity, but the peak (dip) 
density of soliton can be controlled by the Feshbach res- 
onance parameter. The soliton solutions reported here 
may be more realistic and leave scope for more physical 
explanation and application in the future. 



6 



V. ACKNOWLEDGMENTS the NSF of China under Grant No. 10874038, the 

Province Natural Science Foundation of Shanxi under 
This work is supported by Hundred Innovation Tal- Grant No - 2007011007, and the Natural Science Founda- 
cnts Supporting Project of Hebei Province of China, tion of Hebei Province under Grant No. A2008000006. 



[1] N. J. Zabusky, and M. D. Kruskal, Phys. Rev. Lett. 15 
(1965) 240. 

[2] A. Hasegawa, Y. Kodama, Solitons in Optical Commu- 
nications (Oxford University Press, Oxford, 1995). 

[3] M. J. Ablowitz, P. A. Clarkson, Soliton, Nonlinear Evo- 
lution Equations and Inverse Scattering (Cambridge Uni- 
versity Press, Cambridge, 1991). 

[4] F. Dalfovo, S. Giorgini, Lev P. Pitaevskii, and S. 
Stringari, Rev. Mod. Phys. 71 (1999) 463. 

[5] Kevin E. Strecker, Guthrie B. Partridge, Andrew G. Tr- 
uscott & Randall G. Hulet, Nature 417 (2002) 150; 
L. Khaykovich, F. Schreck, G. Ferrari, T. Bourdel, J. 
Cubizolles, L. D. Carr, Y. Castin, C. Salomon, Science 
296 (2002) 1290. 

[6] L. D. Carr and Y. Castin, Phys. Rev. A66 (2002) 063602. 

[7] L. Salasnich, A. Parola, and L. Reatto, Phys. Rev. Lett. 
91 (2003) 080405; 

G. Fibich et al., Phys. Lett. A239 (1998) 167. 
[8] U. Al Khawaja, H. T. C. Stoof, R. G. Hulet, K. E. 
Strecker, and G. B. Partridge, Phys. Rev. Lett. 89 (2002) 
200404. 

[9] P. G. Kevrekidis, G. Theocharis, D. J. Frantzeskakis, and 

Boris A. Malomed, Phys. Rev. Lett. 90 (2003) 230401. 
[10] Lu Li, Zaidong Li, Boris A. Malomed, Dumitru Miha- 

lache, and W. M. Liu, Phys. Rev. A 72 (2005) 033611; 

Lu Li, Boris A. Malomed, Dumitru Mihalache, and W. 

M. Liu, Phys. Rev. E 73 (2006) 066610. 
[11] L. Wu, J.-F. Zhang, and L. Li, New J. Phys. 9 (2007) 69. 
[12] S. L. Cornish, S. T. Thompson, and C. E. Wieman, Phys. 

Rev. Lett. 96 (2006) 170401. 
[13] R. Dum, J. I. Cirac, M. Lewenstein, and P. Zoller, Phys. 

Rev. Lett. 80 (1998) 2972. 
[14] S. Burger, K. Bongs, S. Dettmer, W. Ertmer, K. Sen- 

gstock, A. Sanpera, G. V. Shlyapnikov, and M. Lewen- 
stein, Phys. Rev. Lett. 83 (1999) 5198; 

J. Denschlag et al., Science 287 (2000) 97. 
[15] A. D. Jackson, G. M. Kavoulakis, and C. J. Pethick, 

Phys. Rev. A 58 (1998) 2417. 
[16] Th. Busch and J. R. Anglin. Phys. Rev. Lett. 84 (2000) 

2298; 

C. K. Law, C. M. Chan, P. T. Leung, and M.-C. Chu, 
Phys. Rev. Lett. 85 (2000) 1598; 

Biao Wu, Jie Liu, and Qian Niu, Phys. Rev. Lett. 88 
(2002) 034101. 

[17] J. C. Bronski, L. D. Carr, B. Deconinck, and J. N. Kutz, 
Phys. Rev. Lett. 86 (2001) 1402; 

D. L. Feder, M. S. Pindzola, L. A. Collins, B. I. Schneider, 
and C. W. Clark, Phys. Rev. A 62 (2000) 053606. 

[18] V. N. Serkin, A. Hasegawa, and T. L. Belyaeva, Phys. 

Rev. Lett. 98 (2007) 074102. 
[19] H. H. Chen and C. S. Liu, Phys. Rev. Lett. 37 (1976) 



693; Phys. Fluids 21 (1978) 377 
[20] V. V. Konotop, Phys. Rev. E 47 (1993) 1423; 

V. V. Konotop, O. A. Chubykalo, and L. Vazquez, Phys. 
Rev. E 48 (1993) 563. 

V. V. Konotop, Theoretical and Mathematical Physics, 
99 (1994) 687. 

[21] Qiu-Yan Li, Zai-Dong Li, Shu-Xin Wang, Wei- Wei Song, 
Guangsheng Fu, Opt. Commun. 282, 1676 (2009). 

[22] V. N. Serkin and A. Hasegawa, Phys. Rev. Lett. 85 (2000) 
4502. 

[23] V. N. Serkin, et al., JETP Lett. 74 (2000) 573; 
V. N. Serkin, JETP Lett. 72 (2000) 89; 
V. N. Serkin, et al., IEEE J. Quantum Electron. 8 (2002) 
418. 

[24] Lei Wu, Jie-Fang Zhang, Lu Li, C. Finot, and K. 

Porsezian, Phys. Rev. A 78 (2008) 053807 
[25] J. L. Roberts, N. R. Claussen, James P. Burke, Jr., Chris 

H. Greene, E. A. Cornell, and C. E. Wieman, Phys. Rev. 

Lett. 81 (1998) 5109; 

J. Stenger, S. Inouye, M. R. Andrews, H.-J. Miesner, D. 
M. Stamper-Kurn, and W. Ketterle, Phys. Rev. Lett. 82 
(1999) 2422. 

[26] K. Staliunas, S. Longhi, and G. J. de Valcarcel, Phys. 

Rev. Lett. 89 (2002) 210406. 
[27] P. Engels, C. Atherton, and M. A. Hoefer, Phys. Rev. 

Lett. 98 (2007) 095301. 
[28] F. Kh. Abdullaev, A.M. Kamchatnov, V.V. Konotop, 

and V. A. Brazhnyi, Phys. Rev. Lett. 90 (2003) 230402. 
[29] D. E. Pelinovsky, P. G. Kevrekidis, and D. J. 

Frantzeskakis, Phys. Rev. Lett. 91 (2003) 240201. 
[30] P. G. Kevrekidis and D. J. Frantzeskakis, Mod. Rhys. 

Lett. B 18 (2004) 173. 
[31] V. A. Brazhnyi and V. V. Konotop, Mod. Rhys. Lett. B 

18 (2004) 627. 

[32] Z. X. Liang, Z. D. Zhang, and W. M. Liu, Phys. Rev. 

Lett. 94 (2005) 050402. 
[33] V. M. Perez-Garcfa, V. V. Konotop, and V. A. Brazhnyi, 

Phys. Rev. Lett. 92 (2004) 220403. 
[34] L. Li, Z. H. Li, S. Q. Li, and G. S. Zhou, Opt. Comm. 

234 (2004) 169. 
[35] S. Q. Li, L. Li. Z. H. Li, and G. S. Zhou, J. OPt. Soc. 

Am. B 21 (2004) 2089. 
[36] Z. Y. Xu, L. Li, Z. H. Li, and G. S. Zhou, Phys. Rev. E 

67 (2003) 026603. 
[37] R. Hirota, J. Phys. Soc. Jpn. 51 (1982) 323; 

Zai-Dong Li, Qiu-Yan Li, Xing-Hua Hu, Zhong-Xi Zheng, 

Yubao Sun, Ann. Phys. (N.Y.) 322 (2007) 2545. 
[38] L. Li, Z. H. Li, Z. Y. Xu, G. S. Zhou, and K. H. Spatschek, 

Phys. Rev. E66, 046616 (2002). 



